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SUMMARY 


An inviscid discrete vortex model, with newly derived 
expressions for the tangential velocity imposed at the separation 
points, is used to investigate the symmetric and asymmetric vortex 
separation on cones and tangent ogives* The circumferential 
locations of separation are taken from experimental data. Based on 
a slender body theory, the resulting simultaneous nonlinear 
algebraic equations in a cross-flow plane are solved with Broyden f s 
modified Newton-Raphson method. Total force coefficients are 
obtained through momentum principle with new expressions for 
nonconical flow. It is shown through the method of function 
deflation that multiple solutions exist at large enough angles of 
attack, even with symmetric separation points. These additional 
solutions are asymmetric in vortex separation and produce side force 
coefficients which agree well with data for cones and tangent 
ogives. 
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Lift coefficient 

Lift coefficient at station i 

Normal force coefficient 

Normal force coefficient at station i 

Side force coefficient 

Side force coefficient at station i 

Base area of a body of revolution 

Free-stream velocity 
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Radius of the base area 

Outward normal to a cross section of a body at x 
Nondiraensional lateral coordinate (= y/a) 
Nondimensional vertical coordinate (= z/a) 
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Dimensional vertical coordinate 
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Sideslip angle 
Semi-apex angle of a cone 

Separation angle measured from the starboard 

generator 

Vortex strength 

Nondimensional vortex strength (= r/2iraaU) 
Nondimensional complex coordinate (= Z/a) 
Nondimensional complex conjugate coordinate (■ Z/a) 
Nondimensional complex coordinate of separation 
position 

Nondimensional complex conjugate coordinate of 
separation position 


Velocity potential 



INTRODUCTION 


Recent and new fighters and missiles are developed to have 
increased aerodynamic performance and high maneuverability. These 
requirements frequently involve operation at high angles of attack 
and consequently body separated vortices. Therefore, understanding 
and predicting the flow characteristics near the body can be 
accomplished only if these vortices are accounted for. 

To represent mathematically the body vortex flow field, Bryson 
(Reference 1) introduced two inviscid line vortices in the leeward 
side of the body. He considered symmetric flow past bodies of 
revolution. This approach was extended by Schindel (Reference 2) to 
bodies of elliptic cross section and of cambered longitudinal 
axis. By removing restriction to laterally symmetric flow, Dyer, 
Fiddes, and Smith (Reference 3) extended this simple model to deal 
with asymmetric vortex separation on cones. In all of these 
investigations, a stagnation condition was imposed at the separation 
line as the Kutta condition. On the other hand, Moore (Reference 4) 
imposed a tangential velocity at the separation line based on an 
expression derived by Smith (Reference 5). The effect was found to 
be small in symmetrical flow. 

Other methods used to predict the asymmetric vortex shedding 
include the vortex cloud method (Reference 6) and solving the 
complete Navier-Stokes equations (Reference 7). As described by 
Mendenhall et al. in their paper (Reference 6), in the vortex cloud 
method a perturbation to the symmetric solution must be imposed in 
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order to obtain an asymmetric solution. In fact, every investigator 
has a different means of introducing a perturbation to the symmetric 
solution. The method which Mendenhall et al. used is to modify the 
predicted symmetric separation point on either side of the body so 
that the separation position is not kept symmetric. By solving the 
complete Navier-Stokes equations at high a, the asymmetric vortex 
shedding was found on a cone cylinder (Reference 7). However, the 
asymmetry in the flow field was introduced through asymmetric 
numerical truncation errors with the computed side force being too 
small. 

In the present study, different discrete vortex models will be 
studied in the cross-flow plane. New expressions for the tangential 
velocity at the separation line will be included in the models to 
investigate the vortex flow effect on cones and tangent ogives in 
both symmetric and asymmetric flows. The resulting nonlinear 
algebraic equations are solved for multiple solutions (i.e. both 
symmetric and asymmetric ones) based on a function deflation 
concept. 


MATHEMATICAL MODEL 

Slender-body theory is assumed applicable. This leads to the 
usual representation of the flow in terms of a streamwise 
disturbance velocity which depends only on the streamwise 
coordinate, and a quasi-two-dimensional cross flow which is 
independent of the Mach number. For cones, the flow is assumed to 
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be conical, so that the separation lines are straight lines, OS^ and 
OS 2 (see Figure 1). To represent the nonconical flow, a tangent 
ogive is divided into several stations* Between two stations, flow 
properties are assumed to vary linearly (i.e. 9T/9x = AT/ Ax = 

(I\ - r._ 1 )/Ax, 9a/9x = (a^ - a^-^/Ax, etc.). The separation lines 
are specified by their angular coordinates, 0 ^ and 62 , measured from 
the starboard generator. The sign convention is that a positive 
value of T represents a circulation which is counterclockwise when 
viewed from downstream. With specified separation lines, the 
problem therefore involves six unknown quantities, which include the 
strengths and coordinates of two discrete vortices. These are 
determined by six equations obtained from the following conditions: 
(i) An analogue of the Kutta condition, namely that the 

velocity vectors on the surface of the body at and S 2 
should be directed along the separation lines, OS^ and 
OS 2 • 

In the past investigations, either a stagnation 
condition (Reference 3) or a tangential velocity 
(Reference 4) was imposed at the separation point. Since 
the latter did not provide significant improvement to the 
solution, a new velocity condition is needed. 

In the vortex cloud method (Reference 6 ) the 
relationship between the vortex strength and velocity at 
the boundary layer edge, U^, is used, which is 
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ar/at * k (D 

where K is an empirical factor used to reduce the strength 

of the shedding vorticity* In a subsonic flow, K = 0.6, a 

typical value suggested by most investigators* If half 

value of U is taken as the average tangential velocity at 
e 

the separation point, as assumed in Smith’s model 
(Reference 5), then 


ar 

U — = K £ 2 - 


( 2 ) 


9x 


For a cone with 3T/3x 35 T/x, the tangential velocity 
becomes 

K 2x 2TraxUa v'K 

where y is the nondimens ionali zed vortex strength, and 


y = r/2naaU 


For a tangent ogive, it is assumed that 


8r/3x = Af/Ax, and therefore 
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where i means the present station and i-1 means the 
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previous station. 

However, experimental data (Reference 8) show that 
the vortex convection speed in a shear layer is 
approximately 

V = fU (5) 

tm e 


where f = 0.61, instead of 0.5 used in Equation (2). Then 
for a cone 

V = -^2 Ay<s/a (6) 

tm /K 


and for a tangent ogive 
„ 2fUa 


tm 


/k 


/ir(a i Y i " a i _ 1 Y i _ 1 )/ Axcl 


(7a) 


Six models based on different f and K values are 
investigated for cones in this report. They are 

1. Stagnation separation model (V tm = 0) 

2. Model 1, Smith's separation model (V = Ua/Tryd/ a ) 

3. Model 2, f - 0.61 and K = 1.0 (V = 1.220 Ua/iWa ) 

tm 

4. Model 3, f - 0.61 and K = 0.6 (V = 1.575 Ua/iryd/a ) 

tm 

5. Model 4, f - 1.00 and K = 1.0 (V = 2.000 Ua/Tryfi/a ) 

6. Model 5, f = 1.00 and K = 0.6 (V = 2.582 Ua/iryd/a ) 

tm 

After testing, a best model with f and K values 
written as f^ and will be chosen to investigate the 
asymmetric vortex separation on cones and vortex flow 
separation on tangent ogives. Based on f^ and values, 
V tm ^ or tan S ent ogives becomes 
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(7b) 


V t m 


2f Ua 
b 


/ir(a i Y i - a i _ 1 Y i _ 1 )/Axa 


(ii) A condition based on the fact that the fluid should 
sustain no force* 

In reality, when the flow is separated from a body, 
there will be a vortex sheet with small strength shedding 
from the body and finally rolling up into a concentrated 
vortex. To simplify the physical problem, only the 
concentrated vortex is used here; and a cut which connects 
the vortex core to the separation point is introduced to 
represent the sheet. It is required that the force on the 
cut must be balanced by force on the vortex, as described 
in Bryson* s paper (Reference 1). In the vortex sheet 
model (Reference 9) usually a cut is introduced to connect 
the end of the sheet to the vortex core. The cut is much 


larger in the present model than that in the vortex sheet 
model. Although the force on the cut has always been 
included in the model in the past, a numerical 
experimentation will be conducted in the present study to 
determine its effect on the solution. 


Following the usual development of the slender-body theory, the 
cross-flow velocity potential, <(), is constructed as a solution of 
the two-dimensional form of Laplace’s equation in the plane of 
transverse coodinates y and z (see Figure 1). In a conical flow, 
all such planes are equivalent. The velocity components are the 
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same at points whose y and z coordinates are scaled by x, and the 


vortex circulations are proportional to x. 

The boundary conditions for the far field and on the body are 

( 8 ) 

(9) 


<}> ~ Ua as Z + eo 

z 1 1 


<j> = U6_ 1 

n local 


for I Z I ^iocal X a local 


where U is the free-stream speed, n is the outward normal to the 
cross section of the body at a constant x, <S^ oca -^ is the angle of 
local slope, and &i oca i is the radius of this section. 

The complex potential may be written as 


W 


2 I\ Z - Z 

iUct(Z - ~) + Ua — £nZ + y-r in * 

Z dx Z - (a /Z ) 


Z - Zr 
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Z - (a 2 /Z 2 ) 


( 10 ) 


where the first term comes from the cross-flow stream in the 
presence of a circular body. The second term represents the area 
expansion in the x direction. The third and fourth terms are due to 
the vortices outside the body and their images inside the body. 

Using the nondiraensional variables £ = z/a, y = r/2iraaU, the 
Kutta condition and the force-free condition can be written as 
Kutta Condition 
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where y^ and y 2 are evaluated at the current station i and Ax is a 
dimensional quantity. 

By separating the force-free condition into real and imaginary 
parts, the system of equations which needs to be solved consists of 
six nonlinear algebra equations of six unknowns: y^, y 2 , y ^ , 

z 1 (c i - y l + iZj), y 2 , and z 2 (? 2 = y 2 + iz 2 ). The system can be 
rewritten — with y^, Zj, etc., representing y^_, zj , etc. — as 


FI = y 


2 , 2 . 
+ z i ‘ 1 


1 (cosS^ - y ^) 2 + (sinQ^ - z^ ) 

2 , 2 . 

>2 * ‘2 ~ 1 

If 2 2 2 

(cosS^ - y 2 ) + (sin©^ - z 2 ) 
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( 15 ) 


2 cos 0 . + V /Uct = 0 
1 tra 


F 2 - y 


2 X 2 1 

*1 + Z 1 " 1 


12 2 
(cos02 - y^) + (sinO^ - z^) 


2 X 2 1 

y 2 +z 2 - l 


2 2 2 
(cos 0 2 - y 2 ) + (sin 0 2 - z^) 


2 cos 0 o - V /Ua = 0 
2 tm 
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For a cone, these equations can be simplified to 
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A detailed derivation of these equations is presented in 
Appendices A and B. The last terras in Equations (17)-(20) and (23)- 
(26) are the effect of the force on the vortex cut. These terms are 
eliminated if the force on the cut is not included in the model. 
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During the solution process, FI through F6 may not be zero* As an 
indicator of convergence, a norm based on FI through F6 is defined 
as 

Norm = (FI 2 + F2 2 + F3 2 + F4 2 + F5 2 + F6 2 )^ 2 (27) 

Variables y^ , y^ , y^, z^, y^ and z^ are regarded as a set of 

—8 

solutions when the norm is less than or equal to 10 • 

It is necessary to supplement the Kutta condition as specified 
in (i) above, with the condition that the flow leaves the surface at 
the separation line. This is because the condition that the 
velocity vector lies along a conical ray is satisfied at an 
attachment line as well as at a separation line* The difference 
between the two is that the flow leaves the surface at a separation 
line and approaches it at an attachment line* To confirm the 
solution, the nature of the flow near the separation line is 
checked. When the flow is separated, it would leave the surface to 
downstream side. Therefore the derivative of tangential velocity 
with respect to the angular position 0, dV^_/d0 , should be less than 
zero at a correct separation point. If the derivative is greater 
than zero, the flow is regarded as approaching the body surface 
rather than leaving it. Therefore, the specified separation line is 
a reattachraent line, not a separation line. The expression of the 
derivative is shown in Appendix C. 

After solutions have been found, the local sectional 
coefficients of normal force and side force are calculated by 
integrating the local complex velocity potential in accordance with 
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the momentum principle of fluid mechanics (see Appendix D). The 
total force coefficients can be obtained by integrating the 


sectional coefficients over intervals of longitudinal stations. The 
results are 
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C T = C..cosa 
L N 


( 30 ) 


where ag is the radius of the base area. For a cone, the equations 
for force coefficients can be simplified to 

C = 2«(2Y 1 y 1 U - 2 ‘ " 2 > - 2 V 2 (1 ■ 2 7 2 > + 0 <31) 

+2! y 2 + 2 2 
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A numerical scheme based on Broyden’s concept (Reference 11) is 
used to solve the system of nonlinear algebraic equations. The 
advantages of this scheme are i) unlike in Newton’s method, the 
gradient matrix is calculated by an approximate but inexpensive 
method; for a complicated system of equations, the computing time is 
significantly reduced; and ii) using an optimizer described by 
Broyden, even poor initial guess does not seem to affect the final 
convergence of the method because the optimizer tends to pull the 
values close to the solution. Therefore, the calculation can be 
restarted from where the program stops without convergence. 

In finding the asymmetric solution, a function deflation 
technique is used. The concept of function deflation is described 
as follows (see Reference 12 for details). Let F(X) be the 
nonlinear system of equations to be solved, F representing FI, F2, 
F3, F4, F5, and F6 in the present case and X representing variables 
Y r y 2* V V Z V y 2’ and z 2* Assuming R representing y R > Y R > 
y D 9 Z D > Yd > and z t> is one root of F ( x )> F ( R ) = Dividing 

the original equation F(X) by ||X - R | | , the effect of root R is 
eliminated so that additional roots can be found. Note that | |X - 
R| | represents the norm of X - R. It is defined as 

I |X - K| | - [Tj * r Ri ) 2 + (y 2 - y E2 ) 2 + <y t - y Ei > 2 + 
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(34) 
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Mathematically when X approaches the first root R^ , F(X)/||X - | 

becomes the derivative of F(R^) which does not vanish unless R^ 
happens to be a multiple root. The mathematical expression is as 
follows : 


goo = 

I|X - Rjll 

lira G(X) = lim F ’ (R ) 

X>R 2 X-*R | |X - R | | 


(35) 

(36) 


Experience shows that using the aforementioned concept of 
function deflation, convergence to previously obtained roots can 
always be avoided. 


RESULTS AND DISCUSSIONS 

To investigate the flow on a cone, the separation line must 
first be specified. Since one of the main purposes in the present 
study is to show that the asymmetric vortex separation is part of 
the multiple solutions in the boundary value problem, the symmetric 
separation lines will be based on an empirical equation (Reference 
13) for predicting separation angle in a laminar flow, as 

0 - [73 - (51.4a - 450) 1/2 ] [0.76 + 0.0246] (37) 

Equation (37) has been shown to work very well in the range of 
10° < a < 30° 
for semi cone angles of 
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5° < 6 < 20° 


For tangent ogives, the following equations given in Reference 
14 are used to predict the separation position. That is 

10° < a < 40° 

0 = (13.12 + 6)(3.13 - /0.116a - 1.16) for { (38) 

0° < 6 < 15° 

0° < a < 10° 

0 = 3.13 (13.12 + 6) for { (39) 

0° < 6 < 15° 

As explained by Schindel (Reference 13), the equation of 
predicted separation condition is taken from Crabbe ? s theory which 
agrees with experimental data only in some ranges of angle of attack 
and cone angles* 

In addition, to study the separated flow on tangent ogives, an 
additional information is needed; that is, the axial stations 
downstream of which symmetric and asymmetric vortex separation will 
occur. Numerical experience indicated that at a small a, such as 
three degrees, the flow field can be assumed to be attached and the 
normal force coefficient can be calculated with the slender body 
theory (i.e. C N = 2a). At a slightly higher a, such as four 
degrees, a search is then made to find a station (called x^) at 
which symmetric vortex flow can be calculated: i.e., the boundary 

value problem has a solution. For this purpose, let 

A local = a/6 local (40) 

where <5i oca i is the angle of the local surface slope. The value of 

_ g 

Aiocal at > ^local* * 3e t ^ ie one use ^ at a ii angles of attack 
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to determine where symmetric vortex flow will begin. That is, at 

S 

any station where A^ ocal > ^ocal* s y mmetric vortex flow is assumed 

S 

to exist. For a cone, was determined to be approximately 

1.2. Over the portion of the body with attached flow, the classical 
slender body theory is used to calculate the normal force 
coefficient. 

For asymmetric vortex flow on a tangent ogive, note that 
experimental data indicated that asymmetric flow would not start 
until an a av = 20^ (Reference 15), where 0^ is the cone half angle, 
or, for tangent ogives, 

-i V d 

0 = tan ![ — J ] (41) 

U N /d) - 0.25 

At a = a , a search is made to find the first station after which 
av 

the total side force coefficient (C^.) is nearly zero (a value of 

0.06 is assumed). The value of A i oca ]^ at that station, defined 
a 

as A^ocai* s then used to determine where the asymmetric vortex 
flow will begin at all a > a • For a < & av , only symmetric vortex 
flow is assumed to exist. 

Cones 

Symmetric flow is investigated first, using the six models 
mentioned earlier. Results for symmetric flow are compared with 
experimental data in Figures 2 to 7. The experimental data of 
are taken from Reference 16. In Figure 2 and Figure 3, all flow 
models are with a force on the vortex cut in the force-free 
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condition. In figure 2, the vortex core positions of a five-degree 
cone are presented. It is seen that flow model 4 gives good 
prediction of lateral core positions, but all flow models seem to 
underpredict the vertical positions. Vortex core positions of a 
ten-degree cone are shown in Figure 3. It is found that all models 
overpredict the lateral positions and underpredict the vertical 
positions. This may be explained by a large force on the vortex 
cut. In addition, when the angle of attack is less than 20 degrees, 
the theory does not predict the presence of vortices but experi- 
mental data show otherwise. In Figure 4, results without the force 
on the vortex cut for a five-degree cone are presented. Flow models 
3 and 4 are found to have accurate prediction in both lateral and 
vertical positions. Results for a ten-degree cone for all models 
without the force on the vortex cut in the force-free condition are 
shown in Figure 5. The prediction of vortex core positions is more 
accurate than models with a force on the vortex cut in the force- 
free condition. In all these figures, it is seen that the effect of 
imposing a tangential velocity at separation points is to make core 
positions more inboard and lower as V tm is increased. 

The lift coefficients calculated by all flow models 
with/without a force on the vortex cut in the force-free condition 
are plotted in Figures 6 and 7. Both figures show that is 
reduced by increasing V tm . The force on the vortex cut seems to 
reduce in the stagnation separation model and model 1. However, 
its effect is small in other flow models. Overall, it is found that 
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model 3 without the force on the vortex cut is the best flow model 


which produces results having good agreement with experimental 
data. Therefore, this model will be used to investigate the 
asymmetric vortex separation on cones and the vortex flow on tangent 
ogives. 

After determining a symmetric solution, a function deflation 
technique is used to search for additional solutions. It is found 
that one asymmetric solution at a given a can always be found except 
at low a f s. Pairs of symmetric and asymmetric vortex core positions 
on an eight-degree cone for different incidence parameters A based 
on flow model 3 without a force on the vortex cut in the force-free 
condition are shown in Figure 8. The predicted coefficients of lift 
and side force are compared with experimental data (Reference 16) in 
Figure 9. The prediction of side-force coefficient appears to be 
accurate, but the lift coefficient is overpredicted. This may be 
explained by the inaccuracy of the slender-body theory. It is well 
known that the slender-body theory will overpredict the lift if the 
fineness ratio is not large enough. Therefore, the concept of 
effective angle of attack as shown for delta wings (Reference 17) is 
used to improve the result. When the lift coefficient predicted by 
the slender-body theory is taken to be 2acosa, which is equated to 
the lift coefficient given by a three-dimensional linear theory for 
a body in attached flow, an effective angle of attack, a e , can be 
found. The results predicted by using effective angles of attack 
are shown in Figure 10, The lift coefficient agrees well with 
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experimental data, and the predicted side force appears to be 
slightly lower at high a. The effective angles of attack at 
different geometric a for different bodies are listed in Table 1. 

When the stagnation separation model is used to predict the 
asymmetric separation, the left vortex (with higher vertical 
position) has larger circulation than the right one (with lower 
vertical position). As the tangential velocity, V tm , on the 
separation line is increased, the left vortex will have smaller 
strength than the right one. Therefore, the effect of imposing the 
tangential velocity at the separation points in asymmetric case may 
be regarded as changing the larger vortex strength from one side to 
the other. The vortex strength ratio predicted by the stagnation 
separation model and model 3 are shown in Figure 11. The switch of 
vortex strength can be readily seen in the figure. 

In Moore's paper (Reference 4), a branch of solution was 
regarded as physically meaningless and was called a lower-branch 
solution. Because the vortex core position is lower than the 
separation line, and the derivative of tangential velocity with 
respect to the angular position at the separation line is greater 
than zero, the specified separation line is actually an attachment 
line, not a separation line. Lower-branch solutions are also found 
in asymmetric cases. The pairs of lower-branch symmetric and 
asymmetric vortex core positions are plotted in Figure 12, and their 
strengths are compared with the physically meaningful (upper-branch) 
solutions in Figure 13. The numerical results for the upper and 
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lower branches for an eight-degree cone are listed in Table 2 and 
Table 3. 

Tangent Ogives 

The procedures for calculating vortex flow on tangent ogives 
are described as follows. The symmetric solution at a station at 
which A.i oca i is reached is first determined. This step may require 
several trials to find that particular station. Symmetric solutions 
at stations downstream of this particular station are then 
calculated. These procedures are repeated for all angles of attack. 

To calculate the asymmetric solution, a search is started to 
find the station at which A loca -^ equals A ± oca ^ • Then the converged 
symmetric solution at that station is disturbed to get an initial 
guess of the asymmetric solution* In the present method, the 
magnitudes of y^ and are reduced to make the initial guess 
asymmetric. Through the functional deflation procedure, the 
asymmetric solution can be found. The same idea is applied to all 
downstream stations. 

The converged symmetric solution at a station at which k-^ oca i 
equals A^ ^ is used as the initial guess for the next angle of 
attack. Then all procedures are repeated until all desired angles 
of attack are processed. When a < cx av , only symmetric solution is 
assumed to exist. 

Two tangent ogives with fineness ratios of 5.0 and 3.5 are used 
for correlation. The effect of the number of integration station in 
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the axial direction on the predicted normal force coefficient for a 
tangent ogive of fineness ratio of 5.0 is presented in Figure 14. 

As the number of integration stations is increased, the location 
where the vortex separation first occurs can be more accurately 
predicted. It is seen from Figure 14 that for the normal force 
coefficient, the effect of the number of integration stations is 
small. In addition, the normal force coefficient remains nearly the 
same under the condition of symmetric or asymmetric separation until 
a > 30 degrees. At these higher a f s, the normal force coefficient 
is higher in asymmetric cases. The side force coefficient is shown 
in Figure 15. It is seen that the value of is more sensitive to 

an accurate prediction of where vortex separation first occurs in 
the axial direction. The magnitude of appears to be well 
predicted when the number of integration stations is adequate. At a 
= 30 degrees, the value starts to increase rapidly. However, its 
magnitude is still low and the vortex positions are therefore quite 
close to those in the symmetric case, as shown in Figure 16. In 
Figures 17 and 18 the overpredicted normal force coefficient and 
side force coefficient are improved by the concept of effective 
angle of attack. 

For a tangent ogive of fineness ratio of 3.5, the lift 
coefficient is presented in Figure 19, and the side force 
coefficients in Figure 20. It is seen that the lift coefficient is 
overpredicted. This can again be improved by the concept of 
effective angle of attack, as shown in Figure 21 where the lift 
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coefficient is shown to agree well with experimental data. However, 
the predicted side-force coefficients based on effective angles of 
attack (Figure 22) are small and lower than the experimental data. 
Values of effective angles of attack at different geometric ot’s are 
listed in Table 1. Calculated vortex positions agree well with data 
in Reference 20, as shown in Figure 23. 

Sensitivity of side force coefficients due to asymmetry of 
separation points is illustrated in Figures 24 and 25. In the 
calculation, the left separation point is kept fixed at the value 
used in the symmetric case, and the right separation point is 
shifted up (+A0) and down (-A0) from the symmetric case to obtain 
the asymmetry of separation points. The side force coefficients for 
a cone, being nearly zero, are seen to be insensitive to A0 as shown 
in Figure 24. Since the separation point on the left side is kept 
fixed, the vortex strength is nearly the same as in the symmetric 
case. On the right side, the vortex strength is reduced with the 
separation point shifted up with the resulting vortex being more 
outboard. On the other hand, the side force coefficients for a 
tangent ogive with fineness ratio 5.0 are shown in Figure 25 to grow 
with A0. To produce an experimental side force coefficient of 2.25, 
the right separation point must be shifted up by about 12 degrees. 
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CONCLUSIONS 


Based on modified Bryson's vortex models, symmetric and 
asymmetric vortex separation on cones and tangent ogives were 
investigated. The modifications pertain to a tangential velocity 
imposed at the separation points and a force on the vortex cut in 
the force-free condition. It was concluded that 

1) Adding the force on the vortex cut in the force-free condition 
would push the vortex core more outboard and lower, thus making 
an inaccurate prediction as compared with available data. 

2) Increasing the imposed tangential velocity at the separation 
points would make the calculated vortex core more inboard and 
lower. The strength of the vortex with a higher vertical 
position (i.e. the left one) would decrease, and that of the 
one on the right with lower vertical position would increase in 
asymmetric cases. 

3) Through a numerical scheme based on the concept of function 
deflation, multiple solutions to the problem of vortex 
separation on cones and tangent ogives were shown to exist. 

The asymmetric solutions were shown to exhibit characteristics 
which were observed in experiments with asymmetric vortex 
separation. 

4) The predicted side-force coefficients for cones and tangent 
ogives agreed reasonably well with available data. 

5) Asymmetric separation points could produce vortex solutions 
which generated large enough side forces as measured in 
experiments for tangent ogives, but not for cones. 
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APPENDIX A 


FORMULATION OF KUTTA CONDITION 


The complex velocity potential at a constant x section with 
slope change da/dx can be written as 


2 . r. 

W = -iUa(Z - •§— ) + Ua — Z Z + j-r t 
Z - n Ziri n 

dx 


z - Zj r 2 z - z 2 


2 2iri 2 .-2 

Z - a /Zj Z - a /Z 2 


(Al) 


For a constant slope change in every section (such as a 

cone), da/dx can be replaced by 6. Then 

2 , , T\ 

i + Ua ■ 

Z~ dx 

r. 


dW a x . „ da 1 , 1 r 1 

dZ “ _iUa(1 + U H Z TiiT L Z " - Z. I 2 


— - T ^-1 " 

1 Z - a /Z, 


J 1 


2 r 1 
■nr 




, I, 

2 Z “ a /z 


(A2) 


Using nondimensionalized variables (c and y) such that z - a? 
and T = 2 iraaUy, then 


, , * 2 iraUy, a , 

77 = “iUu(l + —7) + U — — + - — j-t — L - 7 2 , 

dZ ? Z dx ? 2lIi Z Z 1 Z - a /Z 


V-] - 


1 


2 iraaUy 2 


r i \ i = 

\-7 -7 2 - J 


2lli - Z 2 Z - a 2 / Z - 


-iUa(l 


+ \) + U i - iUYjaf— ^ 1 + 

r 2 dx ? 1 ? n «, - 1 


1 ^2 

1 Uy 2 °[t=“c 


2 « 2 - 1 


(A 3 ) 
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JL = -id +1_) +i£ 1_ _ iY r_JL__ _ ?1 i + 

Uadz c 2 dx ?a T^77 j 


’i zz 1 - i 


lY 2 [ Fh 


2 « 2 - 1 


The right-side separation velocity is derived as follows: 

V = i m {-£iK} 

tm 1 a dZ' iOj^ 

Z=ae 

where V tm is the mean tangential velocity and Im is the imaginary 
part of a complex number. Then 


V » Im{-e 1 
tm 1 dZ ’ i0. 


dZ __ i 
Z=ae 


i t/ 

= Im[-e ^of-Kl +ij) + i Y 1 [*jT"r 


C dx 


Z ~ Z 


1 C?i ~ 1 


lY 2 [ r^1 z — ~H] 19 

2 «2-> c-e 1 
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e dx e a 


i6 i 1 c i i6 i 1 

iY l e ^~i0 I i0r 7 1 ~~ iy 2 e [~id~ 


S 6 — 1 


e " C 2 e c 2 - 1 0=0 X 


i0. -i0. . i0 ; r - 1 
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(e - ^Xe 1) 


29 
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10 , 


? 2 ? 2 " 


iy 2 e ‘ 19 ^ i 0 


-}Ua 


(e 1 - ? 2 )(e l l 2 ~ 1) 


= Im{i2cos9 1 " \ ~ + ±y 1 — Jg" 
dx / 1 


' 1 


-i 0 . 


(e - ~ e ) 


iy 


hh " 1 
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= {2cos9^ - 


' 1 
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{2cos9^ - • 
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{2COS0J - Yj 


2 . 2 , 
y l + Z 1 “ 1 


(cos0^ - Yj) 2 + (sinO^ - Zj) 2 


2,2 . 

y 2 + z 2 " 1 


(cos0j - y 2 ) 2 + (sin0^ - z 2 ) 2 


}Uoc 


Define 


XK1 = 


2 . 2 . 
y i + z i " 1 


Ccos6j -y ^ ) + (sinfij - z 


and 


XK2 


2 , 2 , 

y 2 + z 2 " 1 

T 


(cos0^ - y 2 ) + (sin0^ - Z 2 ) 

Then the separation velocity on the right side becomes 
V tm “ U^ 20080 ! ~ Yj * XKl - y 2 * XK2) 


or 


Ua(XKl * y, + XK2 * y 0 - 2cos0.) + V - 0 
1 2 1 tm 

For the separation velocity on the left side. 


" Im { { 


10, 


2 dW 


} 


Z=ae 


Therefore, 


-v tm = {2cos0 2 - y l 


2 . 2 . 
y l + z ! " 1 


(cos 02 - y) 2 + (sin0j - z^) 2 


2 . 2 . 
y 2 + z 2 - 1 
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(A7) 


(A8) 


(A9) 


(A10) 
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Define 


XK3 = 


2 . 2 . 
y l + Z 1 ~ 1 


(cos0 2 - ) 2 + (sin0 2 - 


and 


(All) 


XK4 = - 


2 4 - 2 1 

y 2 + z 2 “ 1 


2 2 

(cos0 2 - y 2 ) + (sin0 2 - z 2 ) 


The separation velocity on the left side becomes 
-V = Ua( 2 c°s0 2 - Y^ * XK3 ~ Y 2 * XK4) 


or 


Ua(-2cos0 o + y, * XK3 + Yo * XK4) - V 

2 l 2 . tm 


= 0 


(A12) 


(A13) 


Define 


FI = XK1 * y, + XK2 * Yo “ 2cos0, + V / cdJ (A14) 

i 2 I tm 

F2 = XK3 * Y, + XK4 * Y 0 " 2cos0„ - V /aU (A15) 

l 2 2 tm 

In Bryson's model, the separation line is required to be a stream 
line of the three dimensional flow, so that V tm , the mean tangential 
cross-flow velocity component at the separation line, is zero. However, 
by considering the behavior of a vortex sheet model at the separation 
line, Smith (Reference 9) derived an expression for V tm , which can be 
applied whether the sheet is present or not. For a body of revolution, 
such as a circular cone, with the separation line along a meridian. 
Smith’s expression becomes 


V = r£ 1 Li 1/2 

tm '-2 - J 
dx 

By the assumption of conical flow, V tm becomes 


(A16) 
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V - r| r ] 1/2 - D[-Xa 2 _]i /2 . u[le ?] 1/2 - WTFJ 7 S <A 17 > 


2Uxair 


2Uaira J 


For nonconical flow, it is assumed that 3r/3x = AT/Ax. Then 

V = [ML] 1/2 = u [ (F i ~ F i~l )a l 1T j 1/2 = u ^ (l ’i ~ r i-l )a j Trq jl/2 
tH! 2 Ax 2UAxa^Tr 2UAxa i 'jra 


- Ua[(a i Y i - * "/Axa ] 1//2 


(A 18 ) 


However, considering the vorticity flux across the boundary layer at 
separation, as in the vortex cloud method (Reference 6), the relationship 
between the vortex strength and the velocity at the boundary layer edge, 

U , may be written as 


3£ _ r ~ 3U e 

3t J 3r Udr 2 


(A 19 ) 


where the nonslip condition is satisfied at r Q . In the vortex cloud 
model, an empirical factor, K is often used to reduce the strength of the 
shedding vorticity. Various investigators recommended a value in the 


range 


of 0.6 < K < 1.0, with 0.6 being a typical value in subsoni< 


flow. It follows that 


1 L = TI 9r = V e 

8t f 


(A 20 ) 


If V tm is equated to one half of the shedding velocity, then for a 


_e = /_U 1 L = JL/U 1 L=M / Tr ' yS/a 
2 2K 3 x /K 2 3 x /K 1fY 


(A 21 ) 


For tangent ogives with 3T/3x - AT/Ax, 
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tm 


/*[ — - — 4 — - — /(a.Y. - a .y j) • n/Axa (A22) 
2K 3x /K 2 9x /K 11 1111 


On the other hand, experimental data (Reference 8) showed that a 
vortex in a shear layer tends to be convected with a speed less than the 
free stream, so that 

V - fU (A23) 

tm e 

where f * 0.61. If this separation condition is used, then V tm for a 
cone becomes 

v . iL /pi . 22° SiWZ (A24) 

c " /f 2 3i /K 


For a tangent ogive, V fcm is 


V 


= 2f /U _ 2fUcx 
tm v'K 2 


(A25) 


9x /Z /( Vi " a i-lVl >Tr/Axa 

Six models are tested for a cone based on these different separation 
conditions. They are 


(1) 

Stagnation 

separation model (V tm = 0) 


(2) 

Model 

1, 

Smith's separation condition 

( V tm = u/lT, Y<5/a) 

(3) 

Model 

2, 

f 

= 0.61, 

and 

K = 1.0 

< V tm * 

1.220 Ua/ifWa) 

(4) 
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= 0.61, 

and 

K = 0.6 

II 

a 

jj 

> 
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. 

II 

and 

K = 1.0 

< v tm ■ 

2.000 Ua/iiWci) 

(6) 

Model 

5, 

f 

o 

o 

• 

•H 

ii 

and 

K = 0.6 


2.585 Ua/SWa) 


After testing, a best model with the best f and K values, f b and K b , 
is chosen to investigate asymmetric vortex separation on cones and vortex 
flow on tangent ogives. Based on f b and K b values, V tm for tangent 
ogives becomes 

2f Ua 

b __ . .. (A26) 
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APPENDIX B 


FORMULATION OF FORCE-FREE CONDITION 

There are two parts of contribution to the total force: 
(i) Due to the inclination of vortex line 


f v - -ipr v [ it. (f - ^ ^L-> - uaz v /dx] 


z+z 


(Bl) 


(ii) Due to potential jump on cut 

F = ipU(dr /dx)(Z - Z ) (B2) 

e v v s 

where Z v is the vortex core position and Z g the separation point 
coordinate. 

The force-free condition is 


F + F = 0, 
v c ’ 


or 


F + F = 0; 
v c 


that is, 
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Dividing both sides by Ua, it is obtained that 
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Since 
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it follows that 
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The force-free condition on the right side (Z y = Z±) becomes 


-i( 1 + — j 
c i 


t. dx 1 - 1 1 


•2 ?1 C 2 - 1 


1 d ?1 dr, 

I — L + I (i - c ) 

a dx ar i dx 1 s i 


(B7) 


and the force-free condition on the left side (Z v = Z 2 ) becomes 
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1 


- 1(1 + + ^ - -^-1 - ^ 


e 2 dx 


C 2 Cj- 1 


^2 - 


1 dc 2 a dr 2 .- - . 

«7^ + ^rr^- (? 2 " « S > 

dx 2 dx 2 


(B8) 


Let 

C 1 = y l + iz l’ c 2 = y 2 + iz 2 * ? s = cos6 l + 1 sin0 i» 
C g = cos0 2 + i sin0 2 . 

Then the force-free condition on the right side becomes 


-i[i 




+ iz^ ) 


7 ] 


1 da 1 y l ~ 12 

° dx y l + i2 l Yl (y l + iz l Hy l ' iz l } " 1 


+ 


r ( y 2 + iz 2 )(y 2 “ iz 2^ “ 1 , 

±Y 2 fyj + iZj “ ~2 ~ iz 2 )[( yi + lz !> (y 2 - iz 2 > - 1] ^ 


iffl 

0 dx 


dr i 

— [ Cy iZj) “ (cos0. - i sin0. ) ] 

1 dx 


0 


-i[l + 


yi - z i 


- i2y lZl 


( y I- 


2 v2 . 2 2 

z x ) + 4y lZl 



da , y l 

- ' ? 2 
dx y x + z x 


2.2 

y, + z i 


) + 



Y l (z l + ly i^ 


(y2+ z 2" 1)[ ^r y 2 > ~ i(z r z 2 )][(y i y 2 + 
2 2 

K y i“ y 2 ) + (zj- z 2 ) ][(y 1 y 2 + z^ 2 - 


vr ~ 1 ( y 2 z r yi z 2 >^ 
1) 2 + (y 2 z 1 - y 1 z 2 )2] 


1 d ^l a d ^l 

o ~ -rs t(y, “ cos0. ) - i(z - sin0 )] = 0 

dx al l dx 1 1 1 1 
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Finally, the above equation can be separated into real and imaginary 


parts to be 


„ y l Z l . 1 da y l , Z 1 

^ , 2 2.2 a - 2 2 Y 1 ' 2 . 2 . 

(y x + z x ) dx Y 1 + z l y i + z l _ 1 

2 2 

Y 2 (y 2 + z 2 “ 0^ z i “ z 2^ y l y 2 + Z 1 Z 2 " ^ + <7i “ y 2 ^" y i z 2 + y 2 z l^ 
U y l - y 2 ) 2 + ( z x " z 2 ) 2 ] [ ^ y i y 2 + z l z 2 ■ 1)2 + (y 2 Z l " y l z 2 )2] 


2 2 

y l ' Z 1 , 1 da Z 1 


. r r 1 . 1 1 1 1 aa 

1 l L 1 + “ 9 TtJ « - 


'12.2 


(y x + z x ) dx y x + z : y l + Z 1 _ 1 


2 , _2 


K y i - y 2 ) + (z 1 - z 2 ) H(y 1 y 2 + z x z 2 - D + (y 2 z x - y 1 z 2 ) ] 


1 d ^1 a dr i 

— = [(y, - COS 0 ) - i(z - sine. )] = 0 

a / aT, 1 1 1 1 

dx 1 dx 


Using a simple difference to represent the derivatives, that is 


da Aa a i a i-l 

I. — - = - — 

dx Ax Ax 


d l _ C 1 ~ c i-l * a i-l /a i 


(y i ' y i-l ' a i-l /a i ) ” i(z i 


z i-l * a i-l /a i } 


dr r i - r i-i 


then the force-free condition on the right side, with yp zj, etc. 
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representing y, , z, , etc. , becomes 
i x i 


-2 


y l Z l 


, 2 . 2.2 a 

(y l + Z l> 


+ 1 S i ~ a i- l y l 


+ Y, 


- 2 2 '1 2 2 

Ax y x + z x y x + Zj ~ 1 


2 . 2 


Y 2^ y 2 + Z 2 " “ z 2 ^ y l y 2 + Z 1 Z 2 “ ^ + ^ y l y 2 ^ _ y l Z 2 + y 2 Z l ^ 

t( y l ~ y 2 ) 2 + (z x - z 2 ) 2 ][( yi y 2 + Zl z 2 - l) 2 + (y 2 - y^) 2 ] 


2 2 
y l ' Z 1 


. a. - a. . z, 

4 r_r, ■ ' i i i 1 i 1-1 1 

M l 9 9 9 J ” « - o o + Y 


2 . 2 n 2 j a 

r l + Z 1 


(y, + Z, ) Ax y7 + z7 * y7 + z7 - 1 


2 2 ' 12.2 

1 + Z 1 y l + Z 1 


2 . 2 


Y 2 (y 2 + z 2 “ 1)[(y i “ y 2^ y l y 2 + Z 1 Z 2 " " (z i “ z 2 )(_y i z 2 + y 2 Z l )] 

Ky, - y 2 ) 2 + («, - z 2 ) 2 H(y 1 y 2 + - l) 2 + (y^ - y^) 2 ] 

a (y l. ‘ y l, . • a i-l /a i ) - 1(2 l * *1. ,• “iW 
i i i-1 l 1^1 

a 

Ax 


Y, - Y, • a. /a. 
i X i X i-1 1 1 1 


ay 


1 . 

i 


Ax 


[(y 1 - cos0^) - i(z 1 - sin6^ )] = 0 


(BIO) 

Similarly the force-free condition on the left side becomes 


-2 


y 2 Z 2 


2 2 S 2 a 

1 2 Z 2 


I a l ~ a i- l y 2 


- Yo 


t u — 2 2 '2 2 2 , 

(y 0 + Z 0 ) Ax y 2 + z 2 y 2 + z 2 “ 1 


2 . 2 


Y 1 ( y l + Z 1 ~ l>[Cz 1 - z 2 )( yi y 2 + Zl z 2 - 1 ) + (y t - y 2 )(-y 1 z 2 + y 2 z l )] 

K y i - y 2 ) 2 + ( z 1 - z 2 ) 2 ][( y 1 y 2 + Zl z 2 - l ) 2 + ( y 2 z A - y ^) 2 ] 


2 2 

y 9 3., "* 3, . Z 0 

‘ H * + - fri ! - - i ^ zl 2 


" Yo 


/ ^ , ^.>. 4 .- — 2 2 *2 2 2 t 

(y 2 + z 2 ^ <*Ax y 2 + Z 2 y 2 + Z 2 " 1 
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Vf! + z i ~ 1)[(y i ~ y 2 )(y i y 2 + z i z 2 " 1} ~ (z i ~ z 2 )( ~ y i z 2 * y 2 z i )] 
[(y l " y 2 )2 + (z l ” z 2 )2jl(y l y 2 + z l z 2 _ 1)2 + (y 2 Z l " y l z 2 )2] 


} - 


(y 2j - y 2, * a i-l /a i } ' 1(Z 2 


Ja “ Za t d . . j c . I , ) 

a i 2 i 2 1-1 1-1 1 2 i 2 i-l 1 ‘ 1 1 


Ax 


a Y, - Y, * a i-l /a i 
i X *1-1 _ 

a 

Y, Ax 
i 


[(y 9 - 


cos0 2 ) - i(z 2 - sin0 2 )] 


= 0 (Bll) 


Define : 


y i z i i a i - Vi 

F3 = -2 — - v -- + 


r l 


, 2 2.2 a 

(y l + z l> 


*" 2 2 

Ax y 1 + z L 


+ Y 


1 


1 2 2 
y l + Z 1 " 1 


Y 2 ( y 2 ^ + z 2 “ " z 2 )(y l y 2 + Z 1 Z 2 " + (y l " y 2 )( “ y l Z 2 + y 2 z l^ 

E ( y i - y 2 ) 2 + - z 2 ) 2 ][(y 1 y 2 + z^ - D 2 + ^ 2 Z \ ~ y i z 2 )2j 


_ y , - y , * a / a Y, " Y, . a /a 

i *i 1 i-l 111 a i A i A i-l i-l' i 


Ax 


Axy, 


(y 1 - cos 0^ ) = 0 


(B12) 


y, " z a - a. 

F4 = -[1 +— i Dl 


, 2 ^ 2.2 J a 

(y i - V 


Ax 


2 . 2 
yi + z i 


+ Y 


'1 


1 2 2 . 
y l + Z 1 " 1 


Y 2 ( y 2 + z 2 DKy! “ y 2 ><y 1 y 2 + z i z 2 ~ " (z i " z 2 )( ~ y i z 2 + y 2 z i^ 

U y l - y 2 ) 2 + (z x - z 2 ) 2 ][(y 1 y 2 + Zl z 2 - l) 2 + (y 2 z x - y^) 2 ] 

z, - z. • a. /a. Y, ~ Y, * ,/a. 

a. 1, 1 1 _ i 1 D D_ I ^ 1 ^ 

— ~ ^ + — — ^ (z. - sin0. ) = 0 

a - a 11 


Ax 


Axy 


1 


( B13) 


F5 = -2 


y 2 Z 2 


, 2 . 2.2 a 

(y 2 + z 2> 


1 a i ~ a l-l y 2 
+ _ o o Y. 


Ax 


y 2 + z 2 


2 2 


y 2 + z 2 ~ 1 
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Y l (y l + z l " 1)I(z i ~ z 2 )(y l y 2 + Z 1 Z 2 ~ + (y i " y 2 ><- y i z 2 + y 2 z l )] 

t( y i ~ y 2 ) 2 + ( Z 1 “ z 2 ^ 2 ^^ y i y 2 + z l z 2 “ + ^ y 2 Z l “ y l Z 2^^ 


• a /a Y, “ Y, • a. ,/a. 

i -i -i-i l -‘ 1 a i *i ‘i-i x -‘ 1 


a. y 2, ' y 2 


a 


Ax 


a 


(y 7 - cos0 9 ) = 0 


Axy, 


(B14) 


y 2 ' Z 2 


1 ^ • ^41 

F6 = -[1 + — f - i -i ■ A- . 1 2 


, 2 x 2 ^ J a * 

( y 2 + z 2 ) Ax 


~ Y 


y 2 + Z 2 


2 2 


y 2 + z 2 - 1 


Y i (y i + z i “ X ^ (y i “ y 2 )(y i y 2 + z^ - O - (z x - z 2 )(- yi z 2 + y 2 z 1 )] 
[( y l - y 2 ) 2 + (z L - z 2 ) 2 ][( yi y 2 + Zl z 2 - l ) 2 + (y 2 z x - y^) 2 ] 


a, z 2. - z 2, * a i-l /a i a, y l - y l. , ' a i-l /a i 

1 i i-1 . i i i-i , . A v 

— + — (z 2 - Sin0 2 ) = 0 


Ax 


Axy, 


(B15) 

For a conical flow, the derivative terms can be approximated by 


da = a ~ ^ _d£ = £ _ dr _ T _ 6f 

- a ’ - ~ x a 

dx x dx x dx 


Therefore, 


F3 = -2 


y l z l 


+ i 


'l 


Y l z l 


f 2 2.2 A 2 2 2 2 

(yj + Zj) y i z i y i z i “ 1 


( y 2 + z 2 - l)[( Zl - z 2 )( yi y 2 + Zl z 2 - 1) + (y x - y 2 )(-y 1 z 2 +y 2 z 1 >] 
t( y i “ y 2 ) 2 + (zj - z 2 ) 2 ] [ ( y^ y 2 + z^ - i) 2 + (y 2 z 1 - y^) 2 ] 

[ y l + (y x - cos0^) ] / A = 0 (B16) 


41 


F4 = -[l + 


y l “ Z 1 
2 2 2 
( y l + Z i> 


] -x 


2 4 . 2 

yi + z . 


Vi 

y. + z, - i 


( y 2 + z 2 “ " y 2 )(y l y 2 + Z 1 Z 2 " " (z l “ z 2 )(_y l Z 2 + y 2 Z l )] 

2 ^ (y l “ y 2 )2 + ^ Z 1 “ z 2^^ y l y 2 + Z 1 Z 2 ” 1 ^ + ^ y 2 Z l ” y l Z 2^^ 

[z^ - (sin0^ - z i ) ] /A = 0 (B 1 7 ) 


F5 = -2 


y? z 


2 2 


, 1 y 2 
+ - -* - Yr 


, 2 ' 2,2 A 2.2 '2 2 2 . 

(y 2 + s 2 > y 2 + z 2 y 2 + z 2 ‘ 1 


(y 2 + zf - DtZj - Z 2 )(y iy2 + z l z 2 ” + ^ y i “ y 2 ^ -y l z 2 + y 2 Z l^ 

1 t(y l “ y 2 )2 + (z l “ z 2 )2j[(y l y 2 + Z 1 Z 2 “ 1)2 + (y 2 Z l ~ y l Z 2 )2j 


[y 2 + (y 2 “ cosQ^JJ/A = 0 


(B 18 ) 


F6 = -[l + 


2 2 
y 2 ' Z 2 


•] -x 


Y 2 y 2 


, 2 . 2 . 2 J A 2 . 2 2 2 . 

(y 2 + z 2 > y 2 + z 2 y 2 + z 2 ‘ 1 


Y i (y i + z i ” 1 ^^ y i " y 2 )( yiy 2 + z i z 2 “ 1 * 
2 2 

t( y l - y 2 ) + (z x ~ z 2 ) IKy^ + z x z 2 

[ z 2 - (sin0 2 - z 2 > ] /A = 0 


(z 1 ~ z 2 )(- yi z 2 + y 2 z 1 >] 

l) 2 + ( y 2 z i - y^) 2 ] 

(B19) 


where A = a/5. 

Define 

Norm = [FI 2 + F2 2 + F3 2 + F4 2 + F5 2 + F6 2 ] 1/2 

Variables y^ y 2> z^ z 2> Yj» and y 2 satisfying the requirement that 

—8 

the Norm be less than or equal to 10 ° are regarded as a set of 
solutions. 

The last terms of Equations (B12) to (B19) are the effect of the 
vortex cut. When testing models without this effect, these terms are 
eliminated. 
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APPENDIX C 


FORMULATION OF THE DERIVATIVE d V tm /d0 


Since the mean tangential velocity on the body is 


V tm ' f 2cos6 ' Y 1 


h'r 1 


+ Y- 


c 2 c 2 - 1 


-ie - } u ° 


f i 0 \ ✓ “”i 0 " v 1 2 / i 0 v f x yj v 

(e - Cj)(e - ? x ) (e - ? 2 )(e - ? 2 ) 


then 


i dv t« . , a , - i>ii<e ' 19 - v« ie - iri ve - ««> 

uS He 2sln9 + Y i 


, i0 N 2 , -10 - .2 

(e - Cj) (e - ? x ) 


/ “ , \ r, 10/ "10 “ x . “i0/ 10 Nn 

(? 2 ?2 " D[ ie (e " ^ 2 ^ “ ie " ^2 ^ 

z i0 v2. -i0 - n 2 

(e - ? 2 ) (e - ? 2 ) 


= -2sin0 + iy^ 


(Cjij - lKCje” 10 - ^e 10 ) 

7 10 v2, -i0 - n 2 

(e - (e - C x ) 


iY 


(c 2 c 2 - l)(c 2 e 16 " ? 2 el0) 
2 TTe 7 T, -ie - .2 

(e - ^ 2 ) (e - ? 2 ) 


(Cl) 


Let 


1 

Ua 


? 1 = Y 1 + ? 1 

? 2 = y 2 + iz 2 * 
It follows that 
dV 

■7 r = _2sine + iY r 
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Kyj+ izj)(yj- izj) - l][(y^+ iz^)(cos9 - i sin9) - (y^- iz^)(cos0 + i sin0)] 
[(cos0 - y^) + i(sin0 - z^)]^[(cos0 - y^) - i(sin0 - z^)]^ 

[(y 2 + ^ z 2^^2 _ ~ Ht(y 9 + iz 9 )(cos0 - isin0) - (y 9 - iz 9 )(cos© + ±sin0) 


2 2 


2 2 ' 


[(cos0 - y 2 ) + i(sin0 - Z 2 )]^[(cos 0 - y^) - i(sin9 - Z 2 )]^ 
2 2 

(y, + z. - l)(z.cos0 - y.sinO) 

= -2sin9 - 2 Yl — ± i ^ + 

[(cos0 - y^ ) + (sin0 - z^) ] 

2 2 

(y 9 + z„ - l)(z 9 cos0 - y 9 sin0) 

2 Y £ —L ■ — - — s— s- (C2) 

[(cos0 - y 2 ) + (sin9 - z 2 ) 1 

A line represents an actual separation line if dV tm /d 0 1 0 < 0 


on the right side and dV tm /d9|g > 0 on the left side. 
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APPENDIX D 


FORMULATION OF THE COEFFICIENTS OF AERODYNAMIC FORCES 

The lateral force is, based on the momentum theory, 

F = / p^(-^ • dA) 

From Reference 10, the local force can be calculated as 


= -2i f dZ + 2i ^ dZ - 2iS Qut (a - ig) + 


W' 


1/2 pU 


out 


in 


2iS . (a - ig) 
m 


(Dl) 


where is the contour of the out-flow plane, with area S Qut , 

D out 

and Ct, is the contour of the in-flow plane, with area S.» . 

“in 

Note that 

F 1 a 2 r 2 a 2 

$ T? dZ = -2ira 2 a rp (Z. - ■—-) + — tt— (Z„ - - ) 

U 1 U I, _ u n 


i Z, 


"i Z, 


1 


2 2 
2ira.r i a 2ita i r 2 a 

= -2Tta^a - 2ira Ua~ (? 1 " + 2wa Ua (? 2 

l i 1 1 ’2 

i i 


) 


= -2ira^a[l + (£j 


i i C 


1 


-) - Y, (C 9 
i i 


1 


■)] 


(D2) 


i i 

To obtain the coefficients of aerodynamic forces, lateral force 
is nondimensionalized by the base area. That is, 

, F . 


C + iC 


(- 


N. 2 ' /9 ..2 
i i;a B 1/2 pU 


) 


(D3) 


For cases in this report, no sideslip effect is assumed. Therefore 
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c + 1C M = ~ 
v. N . 2 

J i 1 ira„ 


{-2Tra^a[l + y l (c^ - ——) - Y 2 (C 


1 1 C, 


'i i t,r 


2ira . all + Y, (?, “ — 

1-1 i-1 i-1 £ 


) - Y- (c 2 ~ r — )] ! - 

, Z i-1 Z i-1 


2i 2 . 21 2 

— "V + ~1 " a i-l a 

Wa B 


= {4a 2 [y, (C, “ l~) " Y 2 ” 

- -• “ i i C, 




U [ Y X (C X - r 1 — ) - r 2 (C 2 - —-) ] + S - 2a^_ 1 } 

■ 1 X i-1 A i-1 C, i-1 i-1 C 0 


C, = y, + z, and C 2 = ?2 + z ? * 

i i i i i i 


y l + iz l 

+ ic -if «y, + it.) - -r— - r 1 ) - 


i a. 


i i i 7, + z, 

i i 


y 2 + iz 


(( + * -t^ii «y, + *«. . 

i i 1 y n + z n 1-1 i-1 i"l 


y 2. + z 2 

1 1 


y, + iz i ^2 + iz 2 

— 2 — 2 ^—} - y 2 {(y 2 + iz 2 ) ” - 2 -i— 2 + 

y? + zf 2 i-l 2 i-l 2 i-l y 2 + z 2 

Vl 1 ±-l i-1 i-1 

2a 2 - 2a 2 _ x } 


— {“ZaJlYj *x (1 2 2 } + Y 2 Z 2 (1 " 2 2 )] + 

i i y, + z, 11 y 0 + z 0 


y l + Z 1 
i i 


y 2 + Z 2 . 
i 1 
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2a i-i ty i V (1 1 7 2 ^ " 

i-l 1-1 y + z. 

i-l i-l 


Y 2 , Z 2 (1 " 2 2 )1} + 

i-l i-l y_ + z, 

i-l i-l 


i-f uajtY, yi a - - y 2 ,y 2 ,d -T T-r>i - 

»B i i + ^ 1 1 \ 


2 a i-l [Y l, /i. , (1 — rr— > - 

i-l i-l y. + z. 

i-l i-l 


Y 2 /2 (1 2 1 2 )] + a i ‘ a i-l } 

i-l i-l y_ + z , 

i-l i-l 


It follows that 


- V J> + 2y 2 *2.° ~ ' 2 ' '• • 2 > + 
1 a B i i y, + z, i i y 2 + z, 

i i i i 


2Yi x z i 1 ~ k r (1 ~ “2 — 7T— > ’ 

i-l i-l a. y, + z. 

i-l i-l 


2 y 2 , , Z 2 . , 1 2 (1 ~1 ~2 

i 1 i 1 3. y 2 z- 

i-1 i-l 


2<Xa i 1 1 

V = ~ {2y l y l (1 "22 } " 2y 2. y 2 (1 2 7 2 } " 

1 a* i 1 y, + z, i i y„ + z„ 


y l + Z 1 
i i 


y 2 + z 2 
i i 


2 V i y i. , -^r (1 5 — TT— > + 

ilila, y. •¥ z , 

i-l i-l 
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2y 


2 y 2 
i-1 i-1 



2 

y 2 


-) + 1 - 


+ z. 


i-1 


'i-i 


‘i-i. 


and 


C 

y 


- I 


i 

= l 


i 



(D7) 


(D8) 

(D9) 


For conical flow (such as on a cone), all flow properties 
changing with respect to x are constant. Therefore, 

C = 2a{-2y 1 z 1 (l — j) + 2y 2 z 2^ 2 ~^ — 2^ (Dl°) 

y y x + y 2 + z 2 

C N = 2a{2y 1 y 1 (l ^ — 2 ) " 2y 2 y 2 (1 — 2 ^ + 1} (D11) 

y l + Z 1 y 2 + Z 2 

C T = C M cosa 
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Table 1: Effective Angles of Attack for Three Bodies 


a = Geometric Angle of Attack 
a g = Effective Angle of Attack 
f.r. = Fineness ratio 



8 Deg. cone 

Tangent Ogive 
of f.r. = 3.5 

Tangent Ogive 
of f.r. = 5.0 

a 

a e 

a e 

a e 

5.0 

— 

4.56 

4.74 

10.0 

8.93 

9.09 

9.43 

12.0 

10.69 

10.88 

11.29 

14.0 

12.44 

12.66 

13.12 

16.0 

14.17 

14.42 

14.95 

18.0 

15.89 

16.16 

16.75 

20.0 

17.58 

17.88 

18.53 

22.0 

19.25 

19.57 

20.27 

24.0 

20.89 

21.24 

21.99 

26.0 

22.51 

22.88 

23.66 

28.0 

24.09 

24.48 

25.31 

30.0 

25.64 

26.05 

26.90 

35.0 

— 

29.78 

30.65 

40.0 

— 

33.16 

— — 


not calculated 
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Table 2: Numerical Results for an Eight-Degree Cone 
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Table 2 (continued) 
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Table 2 (continued) 
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Table 3: Lower-Branch Solutions for an Eight-Degree Cone 
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Table 3 (continued) 
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Root 1.071 0.342 -1.071 


Table 3 (continued) 
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All Models in These Figures Are with f Cut*Terms in Force Free Condition 
Experimental Data (Ref. 13) 

Result from Stagnation Separation Model Result from Model 3 

Result from Model 1 Result from Model 4 

Result from Model 2 Result from Model 5 
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12.0 16.0 20.0 24.0 28.0 32.0 8.0 12.0 16.0 20.0 24.0 28.0 32.0 

c* (Angle of Attack) o< (Angle of Attack) 

Fig.2 Right Vortex Core Positions for Different Models of a 5 Deg. Cone as a 
Function of Angle of Attack (Symmetric Case) 





All Models in These Figures Are with 'Cut* Terms in Force Free Condition 
Experimental Data (Ref. 13) 

Result from Stagnation Separation Model Result from Model 3 

Result from Model 1 Result from Model 4 
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12.0 16.0 20.0 24.0 28.0 32.0 8.0 12.0 16.0 20.0 2 

ot (Angle of Attack) 01 (Angle of Attack) 

Fig.3 Right Vortex Core Positions for Different Models of a 10 Deg. Cone as a 
Function of angle of Attack (Symmetric Cose) 




All Models in These Figures Are without Cut Terms in Force Free Condition 
Experimental Data (Ref. 19) 

Result from Stagnation Separation Model Result from Model 3 

Result from Model 1 Result from Model 4 

Result from Model 2 Result from Model 5 
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12.0 16.0 20.0 24.0 28.0 32.0 8.0 12.0 16.0 20.0 2 

<4 (Angle of Attack) <* (Angle of Attack) 

Fig.4 Right Vortex Core Positions for Different Models of a 5 Deg. Cone as a 
Function of Angle of Attack (Symmetric Case) 




All Models in These Figures Are without Cut' Terms in Force Free Condition 
Experimental Data (Ret. 13) 

Result from Stagnation Separation Modei Result from Model 3 

Result from Model 1 Result from Model 4 

Result from Model 2 Result from Model 5 
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12.0 16.0 20.0 24.0 28.0 32.0 8.0 12.0 16.0 20.0 24.0 28.0 32.0 

<* (Angle of Attack) (Angle of Attack) 

Fig.5 Right Vortex Core Positions for Different Models of a 10 Deg. Cone as a 
Function of Angle of Attack (Symmetric Cose) 





(lift Coefficient) 


All Models i n This Figure Are With ‘Cut’ Terms in Force Free Condition 


G— Experimental Doto (Ref. 18) 

Result from Stagnation Separation Model 

- Result from Model 1 

— — — Result from Model 2 

Result from Model 3 

Result from Model 4 

Result from Model 5 



ot (Angle of Attack) 

FTg.6 Lift Coefficient of o 8 Deg. Cone as a Function of Angle 
of Attack(Symmetric Case) 
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(lift Coefficient) 


All Models in This Figure Are Without 'Cut' Terms in Force Free Condition 


Q- 


Experimentol Data (Ref. 16) 

Result from Stagnation Separation Model 

Result from Model 1 

Result from Model 2 

Result from Model 3 

Result from Model 4 

Result from Model 5 



(Angle of Attack) 

F?g.7 Lift Coefficient of a 8 Deg. Cone as a Function of Angle 
of Attaek(Symmetrlc Case) 
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Result from Model 3 Without 'Cut' Terms In Force Free Condition 



A = 2.500 
Q, - 46.61 


A = 2.750 
6, - 44.66 



A « 2.000 
6, = 51.12 


a 

x a 



FTg.8 Symmetric and Asymmetric Vortex Core Positions 
of a 8 Deg. Cone os a Function of Incidence 
Parameter A (8, +8,- 180 ) 
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Result from Model 3 Without Cut' Terms In Force Free Condition 



FTg.8 Concluded 
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Coefficient) 


G— Experimental Data (Ref. 16) 

Symmetric Result From Model 3 Without 'Cut' Terms 

— — — Asymmetric Result From Model 3 Without 'Cut' Terms 



F1G.9 Uft Coefficient ond Side Force Coefficient of a 8 Deg. 
Cone os a Function of Angle of Attack 
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G— Experimental Data (Ref. 16) 

Symmetric Result From Model 3 Without "Cut" Terms 

- Asymmetric Result From Model 3 Without "Cut" Terms 
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c* (Angle of Attock) 

Fig.10 Uft Coefficient and Side Force Coefficient Calculated at 

Effective Angle of Attack of a 6 Deg. Cone os a Function 
of Angle of Attock 
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I r aa/ r ayl ^° rtex Stren 9 th R<rtio ) ^os/W C Vortex st ™ngth Ratio) 


All Models in These Figures Are Without 'Cut' Terms in 
Force Free Condition 




A (Incidence Parameter) 

Fig. 11 Vortex Strength Ratio for Different Models of a 
8 Deg. Cone as a Function of Incidence 
Parameter A 
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Result from Model 3 Without Cut' Terms in Force Free Condition 


□ 



Fig. 12 Low Branch Symmetric and Asymmetric Vortex Core 
Positions of a 8 Deg. Cone as a Function of 
Incidence Parameter A (0, +0 4 » 180 ) 
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Result from Model 3 Without Cut' Terms in Force Free Condition 


s 


B 



A - 3.000 
6, = 42.85 



A - 3.250 
0, = 41.15 


Fij.12 Concluded 
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Model 3 Without Cut 'Terme In Force Free Condition 

Up Branch Reeult 

— — — Low Branch Reeult 




FTg.13 Vortex Strength Ratio for Both Branchs of a 8 Deg. 
Cone as a Function of Incidence Parameter A 


70 





(Normal Force Coefficient) 


6 

©— Experimental Doto (Re -0.8 *10 .M-0.70, Ref.18) 

------ Result from Model 3 for 25 Stations (Symmetric Cose) 

Result from Model 3 for 25 Stations (Asymmetric Cose) 



Flg.14 Normal Force Coefficient of o Tangent Ogive (Fineness Ratio - 5.) 
as o Function of Angle of Attack 
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(Normal fi 


©— Experimental Data (Re “0.8 *10 ,M*>0.70, Ref. 18) 

Result from Model 3 for 50 Stations (Symmetric Cose) 

Result from Model 3 for 50 Stations (Asymmetric Cose) 



0.0 5.0 10.0 15.0 20.0 25.0 30.0 35.0 40.0 

ot (Angle of Attack) 

F!g.14 Concluded! 
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|Cy| (Side Force Coefficient) 


G— Experimental Data (Re "0.8 • 10 ,M»0.70, Ref.18) 

Result from Model 3 for 50 Stations 

Result from Model 3 for 25 Stations 



Flg.15 Side Force Coefficient of a Tangent Ogive (Fineness Ratio - 5.) 
os a Function of Angle of Attack (Asymmetric Case) 
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Result from Model 3, For 25 Stations (Symmetric Case) 
Result from Model 3, For 25 Stations (Asymmetric Case) 




Fig. 16 Tracing of Vortex Core Position In o Tangent Ogive 
(F in e n ess Ratio « 5). 
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O— Experimental Data (Re -0.8 *10 ,M=0.70, Ref.18) 

------ Result from Model 3 for 25 Stotions (Symmetric Cose) 

Result from Model 3 for 25 Stations (Asymmetric Cose) 



FTg.17 Normal Force Coefficient Calculated at Effective Angle of 
Attack of a Tangent Ogive (Fineness Ratio « 5.) as a 
Function of Angle of Attack 
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(Normal Force Coefficient) 


G— Experimental Data (Re “0.8 *10 ,M*=0.7Q, Ref. 18) 

Result from Model 3 for 50 Stations (Symmetric Cose) 

Result from Model 3 for 50 Stations (Asymmetric Case) 



£* (Angle of Attock) 
FTg-17 Concluded 
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|CyJ (Side Force Coefficient) 


O— Experimental Goto (Ft *0.8 • 10 ,M=0.70, Ref. 18) 

Result from Model 3 for 50 Stations 

Result from Model 3 for 25 Stations 



Fig. 18 Side Force Coefficient Calculated at Effective Angle of 
Attack of a Tangent Ogive (Fineness Ratio * 5.) as a 
Function of Angle of Attack (Asymmetric Case) 
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0— Experimental Data (Re • 0.35 *10 Ref.16) 

t>— Experimental Data (Re “O.SO *10 6 > M**0.80, Ref.19) 

Result from Model 3 for 35 Stations (Symmetric Cose) 

Result from Model 3 for 35 Stations (Asymmetric Cose) 



Fig. 19 Normal Force Coefficient of a Tangent Ogive (Fineness Ratio 3.5) 
as a Function of Angle of Attack 
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(Norma! Force Coefficient) 


O— Experimental Doto (R* =0.35 *10* Ref.16) 

>— Experimental Doto (Re =0.80 •10*M»0.8 I Ref.19) 

Result from Model 3 for 50 Stations (Symmetric Cose) 

Result from Model 3 for 50 StotJons (Asymmetric Case) 



FTg.19 Concluded 
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G- 

>- 


Experimental Data (Re =0.35 • 10 , Ref.16) 

Experimental Data (Re =0.80 • 10 6 ,M=0.8, Ref. 19) 
- Result from Model 3 for 50 Stations 
* Result from Model 3 for 35 Stations 



FTg.20 Side Force Coefficient of o Tongent Ogive (Fineness Ratio = 3.5) 
as a Function of Angle of Attock (Asymmetric Case) 
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(Normal Force Coefficient) 


Q- 

>- 


Experimental Data (Re *0.35 *1Q 6 Ref.16) 

Experimental Data (Re -0.80 *10 6 > M»0.8, Ref. 19) 

- Result from Model 3 for 35 Stations (Symmetric Cose) 

- Result from Model 3 for 35 Stotions (Asymmetric Cose) 



Flg.21 Normal Force Coefficient Calculated at Effective Angle of 
Attack of a Tangent Ogive (Fineness Ratio - 3.5) as a 
Function of Angle of Attack 
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(Normal Force Coefficient) 


i 


G— Experimental Data (Re =0.35 *10 Ref.16) 

>— Experimental Data (Re =0.80 *10tM=0.8. Ref.19) 

Reeult from Model 3 for 50 Stations (Symmetric Cose) 

Result from Model 3 for 50 Stations (Asymmetric Case) 



F1g.21 Concluded 
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|Cy| (Side Force Coefficient) 


G— Experimental Data (R* *0.35 • 10 g , Rtf.16) 

t>— Experimental Data (fe -0.80 • 10,M=0.8, Ref. 19) 

Result from Model 3 for 50 Stations 

------ Result from Model 3 for 35 Stations 



F1g.22 Side Force Coefficient Calculated at Effective Angle of 
Attack of a Tangent Ogive (Fineness Ratio - 3.5) as a 
as a Function of Angle of Attack (Asymmetric Case) 
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- 10.0 - 8.0 - 6.0 



A © (-© . - 6 ) 

° 1 right toft ' 

Fig. 24 The Sensitivity of Vortex Strength and Side Force 
Coefficient Due to The Asymmetry of Separation 
Points on an 8 Degree Cone 
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-10.0 -8.0 -6.0 -4.0 -2.0 0.0 2.0 4.0 6.0 8.0 10.0 

A9 <-V*-«W 


Ftg. 25 The Sensitivity of Side Force Coefficient Due to 

The Asymmetry of Separation Points on o Tangent 
Ogive (Fineness Ratio - 5.0) 
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principle with new expressions for nonconi cal flow. It is shown through the 
method of function deflation that multiple solutions exist at large enough 
angles of attack, even with symmetric separation points. These additional 
solutions are asymmetric in vortex separation and produce side force 
coefficients which agree well with data for cones and tangent ogives. 
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